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Abstract 

Background/Objectives: Recent large-scale genome-wide association studies have identified multiple loci robustly 
associated with BIVll, predominantly in European ancestry (EA) populations. However, associations of these loci with obesity 
and related traits have not been well described in Chinese Hans. This study aimed to investigate whether BMI-associated 
loci are, individually and collectively, associated with adiposity-related traits and obesity in Chinese Hans and whether these 
associations are modified by physical activity (PA). 

Subjects/Methods:\Ne genotyped 28 BMI-associated single nucleotide polymorphisms (SNPs) in a population-based cohort 
including 2,894 unrelated Han Chinese. Genetic risk score (GRS), EA and East Asian ancestry (EAA) GRSs were calculated by 
adding BMI-increasing alleles based on all, EA and EAA identified SNPs, respectively. Interactions of GRS and PA were 
examined by including the interaction-term in the regression model. 

Results: Individually, 26 of 28 SNPs showed directionally consistent effects on BMI, and associations of four loci {TMEM18, 
PCSKl, BDNF and MAP2K5) reached nominal significance (P<0.05). The GRS was associated with increased BMI, trunk fat and 
body fat percentages; and increased risk of obesity and overweight (all P<0.05). Effect sizes (0.11 vs. 0.17 kg/m^) and 
explained variance (0.90% vs. 1.45%) of GRS for BMI tended to be lower in Chinese Hans than in Europeans. The EA GRS and 
EAA GRS were associated with 0.1 1 and 0.13 kg/m^ higher BMI, respectively. In addition, we found that PA attenuated the 

effect of the GRS on BMI (Pir,teraction = 0.022). 

Conclusions: Our observations suggest that the combined effect of obesity-susceptibility loci on BMI tended to be lower in 
Han Chinese than in EA. The overall, EA and EAA GRSs exert similar effects on adiposity traits. Genetic predisposition to 
increased BMI is attenuated by PA in this population of Han Chinese. 
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Introduction 

Obesity is a major rislc factor for chronic diseases including type 
2 diabetes, hypertension, and cardiovascular diseases [1], and 
results from the interaction between genes and environmental 
factors [2] . Recent genome-wide association studies (GWAS) have 
identified more than 50 loci that are convincingly associated with 
BMI in populations of European ancestry [3-7] and four 
additional loci {PCSKl, CDKALl, KLF9 and GP2) in East Asian 
populations [8-9] . However, Asians in general have a higher body 
fat percentage than Europeans at any given BMI [10], and the 
transferability of loci identified in European ancestry populations 



to Chinese populations has not been extensively studied. Although 
1 8 loci identified in European populations (FTO, SECl 6B, MC4R, 
GIPR, RBJ, BDNF, MAP2K5, GNPDA2, TFAP2B, TMEM18, 
TMEM160, TNNBK, MTCH2, FAIM2, FLJ35779, SH2B1, 
RPL27A and MEGRT) have been confirmed in East Asian 
populations [9], only 12 BMI-associated single nucleotide poly- 
morphisms (SNPs) identified before 2010 have been investigated 
for their effects on BMI in Chinese populations [11-16], 
Furthermore, the collective contribution of these genetic loci to 
adiposity-related traits has not been studied in Chinese. Therefore, 
we aimed to examine whether individual and combined SNP 
effects, evaluated by genetic risk score (GRS), are associated with 
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BMI and adiposity-related traits in Cliinese Hans, and whether 
there are difierences in genetic predisposition between Chinese 
Hans and European ancestry populations. Previous studies have 

shown that increased physical activity (PA) can attenuate the 
genetic predisposition to increased BMI and obesity risk [17—26]. 
Thcrcfor(;, we also examined whether there were interactions 
between PA and GRS on BMI in Chinese Hans. 

Materials and Methods 

Study Population 

The Nutrition and Health of Aging Population in China 
(NHAPC) study is a population-based study of 3,210 unrelated 
Chinese Hans, aged 50-70 years living in Beijing or Shanghai. A 
total of 2,894 individuals had DNA samples are available for 
genotyping. The study have been described in detail before [27]. A 
complete physical examination was conducted in local health 
stations or community clinics, during which body weight and 
height were measured, and overnight fasting blood sample were 
collected. PA was assessed by a standardized questionnaire 
(International Physical Activity Questionnaire, short last 7-day 
format). Physical activity level for each individual was classified as 
low, moderate, and high according to the questionnaire protocol 
(http://www.ipaq.ki.se/scoring.pdf). BMI was calculated as 
weight/height^ (kg/m^). Normal weight, overweight and obesity 
were defined as BMK24, 24<BMI<28 and BMI>28 kg/m^ 
respectively, using the BMI cut-off points for Chinese populations 
[28]. In addition, 1,016 Shanghai participants completed whole- 
body DXA scan by Hologic QDR 4500 W scanner (Hologic, 
Bedford, MA). Total l)ody fat, trunk fat and leg fat, derived using 
the software provided by the manufacturer, were divided by total 
body mass to generate their respective percentages in whole body. 
The phenotypic characteristics of the population are shown in 
Table SI. The study protocol was approved by the Institutional 
Review Board of the Shanghai Institute for Nutritional Sciences, 
and all participants provided written informed consents. 

SNP Selection 

Of more than 60 obesity-related SNPs established to date, we 
focused on the 36 loci identified in GWAS for BMI in European 
and East Asian ancestry populations [3-9]. We selected one 
representative SNP for each locus. Since eight SNPs were 
monomorphic in Chinese Hans, a total of 28 BMI-associated 
SNPs were included in the current analyses, of which 24 SNPs 
were first identified in European ancestry populations and four 
SNPs were first identified in East Asian populations. 

Genotyping 

Genotyping and imputation of SNPs have been described in 
detail before [29]. In brief, all samples were genotyped using the 
lUumina Human660W-Quad BeadChip (lUumina, Inc., San 
Diego, CA) and sample call nitcs exceeded 97%. Each SNP had 
a call rate >95"/() and genotype distribution was in Hardy- 
Weinberg erjuilibrium P>0.01. We used the IMPUTE software to 
impute SNPs from the phase 2 HapMap CHB+JPT (release 
number 22) reference panel (http://mathgen.stats.ox.ac.uk/ 
impute/impute.html, version 2.1.2), and the 11 imputed SNPs 
had a high quality (Proper_info>0.8) (Table S2). 

Genotype Risk Score Calculation 

Genotypes were coded as 0, 1 , or 2 according to the number of 
BMI-increasing alleles of each SNP, and BMI-increasing alleles 
were defined according to findings of previous GWAS [3-9] . The 
overall GRS was calculated by summing BMI-increasing alleles of 



all 28 SNPs. In addition, we calculated a EA GRS by adding the 
BMI-increasing alleles of the 24 SNPs that reached genome-wide 
significance (P<5xl0 ^) in European ancestry populations [7], 
and a EAA GRS by adding the 1 1 BMI-associated SNPs that 
reached genome-wide significance reported by Wen et al and 
Okada et al [8-9]. Seven SNPs are shared by EA and EAA GRS 
(Table S2). Missing genotypes were replaced by the average allele 
count of respective SNP to calculate GRS. No individuals had 
more than three missing genotypes. 

Statistical Analysis 

We used general linear regression models to test the associations 
of individual BMI-related SNPs (both the genotyped and imputed) 
and GRSs with adiposity-related traits, including BMI, body fat 
percentage, trunk fat percentage and leg fat percentage, assuming 
an additive effect of the BMI-increasing alleles. In addition, we 
inverse-normally transformed all outcome variables (to a mean of 
0, and SD of 1) and performed the same analyses. This allows us to 
compare per- allele effect sizes of individual SNPs and of the GRSs 
across BMI, body fat percentage, trunk fat percentage and leg fat 
percentage. Logistic regression models were used to test associa- 
tions between individual SNPs or GRSs and risks of obesity and 
overweight. Genotype distributions of European ancestry and 
Chinese Hans were compared using a likelihood ratio test. 
Heterogeneity of effect sizes between Chinese and Europeans 
was calculated with Cochran's Q;test. The explained variance for 
BMI was calculated as 2*f*(l-f)*fi^, where f is the minor allele 
frequency and (3 is the per-allele effect on standardized \ alu(;s of 
BMI. Interactions between GRS and PA were tested by including 
GRS, PA, and the interaction terms (GRS*PA) in general linear 
regression model. The associations of genetic variants with BMI 
and obesity and interaction analyses were adjusted for age, age'^, 
sex, region (Beijing and Shanghai), and the first two principle 
components. The associations of genetic variants with body, trunk 
and leg fat percentage were adjusted for age, age^, and sex. The 
power calculations for association analyses were performed using 
Quanto software (available at http://hydra.usc.edu/gxe/) using 
the minor allele frequencies and mean values of BMI from the 
NHAPC study and the effect sizes originally reported by previous 
studies [3-9] . Tlu- power calculations for interaction analyses were 
performed using code-generating program MLPowSim to gener- 
ate R code for simulations in interaction analyses [30]. AH tests 
were two-sided, with nominal significance defined as /^0.05. 
Bonferroni correction was used to adjust for multiple testing in the 
()l)esity related traits analyses, and -PsO.0004 was considered to be 
statistical significance. Statistical analyses were performed in R 
software (version 2.13.1; www.r-project.org). 

Results 

Associations with BMI for 26 of the 28 SNPs (binomial test: 
P= 1 .52 X 10 ) were directionally consistent with those previously 
reported in Europeans and East Asians, of which four, in or near 
TMEM18, PCSKl, BDNF and MAP2K5, reached nominaUy 
significance (P<0.05) (Table 1). Two loci, TFAP2B identified in 
Europeans and KLF9 identified in East Asians only, showed (non- 
significant) opposite effects on BMI compared to those reported 
previously. The BMI-increasing alleles were nominally associated 
with higher body fat percentage for one SNP in MTIF3 (binomial 
P— 0.002), with greater trunk fat percentage for two SNPs in/near 
MTIFS and SH2B1 (binomial /"= 0.006) and with less leg fat 
percentage for one SNP near PCSKl in 28 SNPs (binomial 
P= 0.092) (Table 2). Three BMI-increasing alleles of SNPs in or 
near TMEM18 and PCSKl were nominally associated with 
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Table 1. 


Associations of individual SNPs and GRSs with BMI. 










Gene 


SNP 


EAF 


BMI (kg/m^) 




Explain variance (%) 


Power 








Beta (SE) 


P 






NEGRI 


rs2568958 


0.91 


0.02 (0.16) 


0.90 


0.00 


0.20 


TNNI3K 


rsl514175 


0.78 


0.10 (0.11) 


0.38 


0.03 


0.15 


PTBP2 


rsl 555543 


0.88 


0.14 (0.14) 


0.30 


0.02 


0.11 


SEC168 


rs574367 


0.20 


0.19 (0.12) 


0.11 


0.12 


0.59 


TMEM18 


rsl 11 27485 


0.91 


0.43 (0.16) 


0.0083 


0.20 


0.59 


RBJ 


rs6545814 


0.42 


0.10 (0.09) 


0.30 


0.03 


0.43 


ETV5 


rs7647305 


0.95 


0.35 (0.20) 


0.087 


0.09 


0.16 


GNPDA2 


rsl 0938397 


0.30 


0.10 (0.10) 


0.33 


0.03 


0.54 


FU3S779 


rs21 12347 


0.43 


0.06 (0.10) 


0.52 


0.01 


0.27 


PCSKl 


rs261967 


0.42 


0.23 (0.09) 


0.014 


0.21 


0.81 


CDKALl 


rs9356744 


0.40 


0.13 (0.09) 


0.16 


0.07 


0.38 


NUDT3 


rs206936 


0.53 


0.12 (0.09) 


0.18 


0.04 


0.16 


TFAP2B 


rs987237 


0.16 


-0.14 (0.12) 


0.28 


0.04 


0.34 


LRRN6C 


rsl 0968576 


0.22 


0.11 (0.11) 


0.31 


0.05 


0.25 


KLF9 


rsl 1142387 


0.33 


-0.005 (0.10) 


0.96 


0.00 


0.09 


RPL27A 


rs4929949 


0.41 


0.08 (0.09) 


0.37 


0.02 


0.15 


BDNF 


rsl0501087 


0.53 


0.18 (0.09) 


0.050 


0.15 


0.64 


MTCH2 


rs3817334 


0.32 


0.07 (0.10) 


0.46 


0.03 


0.15 


FAIM2 


rs71 38803 


0.28 


0.17 (0.10) 


0.10 


0.09 


0.31 


MTIF3 


rs4771122 


0.19 


0.23 (0.13) 


0.067 


0.09 


0.18 


MAP2K5 


rs4776970 


0.24 


0.21 (0.11) 


0.048 


0.12 


0.32 


GP2 


rsl 2597579 


0.29 


0.04 (0.10) 


0.67 


0.00 


0.10 


SH2B1 


rs4788102 


0.14 


0.07 (0.13) 


0.61 


0.00 


0.29 


FTO 


rs9939609 


0.11 


0.12 (0.14) 


0.40 


0.03 


0.83 


MC4R 


rsl 778231 3 


0.23 


0.01 (0.11) 


0.95 


0.00 


0.66 


KCTD15 


rs29941 


0.24 


0.05 (0.11) 


0.65 


0.01 


0.13 


GIPR 


rsl1671664 


0.54 


0.04 (0.09) 


0.67 


0.00 


0.47 


T MEM 160 


rs3810291 


0.29 


0.19 (0.10) 


0.074 


0.12 


0.22 


GR5 




23.4±3.2 


0.11 (0.02) 


1 .54E-07 


0.90 


0.99 


EA GRS 




19.2±2.9 


0.11 (0.02) 


2.11E-06 


0.76 


0.99 


EAA GRS 




8.4±2.1 


0.13 (0.03) 


6.79E-05 


0.56 


0.99 



Data are beta (SE) per BMI-increasing allele, 
doi:1 0.1 371 /journal.pone.0091 442.t001 



adjusted for age, age^, sex, region and the first two principle components. EAF: effect allele frequency. 



increased risk of obesity (OR was 1.45 and 1.23, /"= 0.002 on a 
binomial test) (Table 3). No associations of individual SNPs with 
adiposity related traits or obesity reached statistical significance (all 
Pvalues>0.004). 

Compared to the effect sizes for the BMI-loci reported in 
European-ancestry populations by Spehotes et al [7] , we observed 
significant inter-ancestry heterogeneity for the effect of TFAP2B 
(■Pheterogeneity<0.05), witli thc cflfcct beuig larger in Europeans than 
in Chinese (Table 1, Figure lA). The BMI-increasing allele 
frequencies of six SNPs were significantly higher and those of 
another 19 SNPs were lower in Chinese Hans than in European 
ancestry populations {P for heterogeneitySO.006), whereas for 
three SNPs, effect allele frequencies were similar (Table 1, 
Figure IB). The explained variances of the four loci that showed 
nominally significant association with BMI {TMEAI18, PCSKl, 
BDMF and MAP2K5) were larger (0.12-0.21% versus 0.01-0.15%) 
in Chinese than in Europeans (Table 1, Figure IC). 



The GRS, which assesses the combined effect of the 28 SNPs, 
showed significant association with adiposity-related traits. Each 
additional BMI-increasing allele in the GRS was associated with 
0.11 kg/m-' higher BMI (P= 1.54x10"'; equivalent to 318 g/ 
allele for a person of 1.7 m) (Table 1), 0.14% higher body fat 
percentage (/"= 0.0090), and 0.10% higher trunk fat percentage 
(P= 0.0082), but not with leg fat percentage (P=0.15) (Table 2). 
The effect and explained variance of GRS on BMI tended to be 
smaller in Chinese (0.11 kg/m^; 0.90%)) than that in European 
ancestry (0.17 kg/m^; 1.45%). 

Using inverse normally transformed BMI-related traits, we 
found that the GRS of 28 SNPs showed the strongest association 
with BMI, followed by trunk fat and overall body fat percentage, 
whereas no association was observed with leg fat percentage 
(Table S3). The GRS was associated with increased risk of obesity 
(OR 1.06, 95% CI 1.03-1.10; P= 0.0006), and overweight (OR 
1.04, 95% CI 1.01-1.06; P= 0.0017) (Table 3). 
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Figure 1. Correlation of BMI-related loci between European ancestry and Chinese Hans. The plots showed the effect size (kg/m^) of loci 
on BMI (A), BMI-increasing allele frequencies (B), and explained variance (%) of loci on BMI (C) in GIANT (Speliotes et al) and NHAPC (our data) studies. 
X-axis represents the data from our study (NHAPC) and Y-axis represents the data from European ancestry populations (GIANT). 
doi:1 0.1 371 /journal.pone.0091 442.g001 



Each additional BMI-increasing allele in the European ancestry 
(EA) GRS, of 24 SNPs, was significantly associated with 0. 1 1 kg/ 
higher BMI {P= 2.11x10 whereas the East Asian ancestry 
(EAA) GRS, of 1 1 SNPs, was associated with 0.13 kg/m^ higher 
BMI (P=6.79xl0""^). In addition, both of EA GRS and EAA 
GRS were associated with increased risk of obesity (OR [95% CI]: 
1.06[1.02-1.10]; P= 0.0044 and 1.10[1.04-1.16]; P= 0.0009, 
respectively) and overweight (OR [95% CI]: 1.03 (1.01-1.06); 
P= 0.015 and 1.06 (1.02-1.09); P= 0.0038, respectively) (Table 3). 
The EA GRS was significantly associated with body fat 



distribution including total fat, trunk fat, and leg fat percentage 
while the associations for EAA GRS did not reach significance 
(Table 2, Figure 2). However, no heterogeneity was found between 
the effects of EA and EAA GRS on adiposity traits {P for 
heterogeneity>0.05) (Figure 2). 

Increased PA levels significantly attenuated genetic effect of 



GRS on BMI (beta; 



■interaction ~ 



-0.06 kg/m^, Pi 



interaction ' 



= 0.022; 



of the GRS on BMI between different levels of PA; Table 4). More 
specifically, the effect of the GRS on BMI was 60% smaller in 
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♦ 

BMI ^ 



Body fat percentage 

♦ 
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♦ 

Trunk fat percentage ; 



. ♦ 

Leg fat percentage ^ 


' ' ' I I I ■ ■ 

-0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05 

Figure 2. The associations of overall GRS, EA GRS and EAA GRS 
with adiposity traits. Figure shows the effects of overall GRS (top 
estimate), EA GRS (middle estimate), and EAA GRS (bottom estimate) on 
inverse normal transformed BMI, body fat percentage, trunk fat 
percentage, and leg fat percentage in our population. 
doi:1 0.1 371 /journal.pone.0091 442.g002 

people with high physical activity (beta = 0.08 kg/m^ per aUele; 
P= 0.004) compared to those with low PA (beta = 0.20 kg/ per 
allele; P= 0.009) (Table 4). The interaction effects were similar for 
the EA GRS (betai„,,„,t;„„ = -0.05 kg/m^; P=0.13) and EAA 
GRS (betainteraction = "O.OS kg/m^; P= 0.044), but did not reach 
nominal significance in the former. Individually, TMENIW- 
rsl 1127485 was nominally significant interacted with PA on 
BMI (beta;„,„„,i„n = 0.53 kg/m^; P=0.04) (Table S4). After 
SEC] 6B-rs57 4367 and FrO-rs9939609 which were reported 
interacted with PA on BMI previously [26], and TMEM18- 
rsl 1127485 were excluded from the GRS, the interaction of 
combined genetic variants with PA on BMI remained significant 
(beta;„,„,,,i<,„= -0.07 kg/m^; P=0.04). 

Discussion 

In this population-based cohort study of Han Chinese, the 
majority of the previously identified BMI loci showed directionaUy 
consistent effects on BMI, of which four reached nominal 
significance. The GRS, based on all 28 SNP, was significandy 
associated with BMI, body fat percentage and trunk fat percentage 
as well as increased risks of obesity and overweight. Effect size and 
explained variance of GRS for BMI tended to be lower in Chinese 
Hans than in European ancestry populations [7]. The overall 
GRS, EA GRS and EAA GRS had similar effects on adiposity 
traits. Furthermore, PA attenuated the association between GRS 
and adiposity related traits; the effect of GRS on BMI was 
significantly larger in individuals with low compared to those with 
high physical activity level. 

We observed that associations of most loci studied in our Han 
Chinese population were directionaUy consistent with those 
reported for European and East Asian ancestry populations [3- 
9] . Nevertheless, the relative impact of the EA-identified loci in 



our population tends to be smaller than in European ancestry 
populations. For example, the FTO SNP, which has the largest 
effect size (0.39 kg/m^ per allele) and explained variance (0.34%) 
for BMI in European ancestry populations [5-7], has only a 
modest effect (0.12 kg/m^ per aUele; 0.03%) on BMI m our 
population. Previous finding from meta-analysis conducted in 
~ 10,000 East and South Asians also provided evidence that the 
effect of FTO on BMI is smaller than that observed in European 
ancestry [31]. We also found nominally significant ancestry 
differences in TFAP2B that effect size was smaller in Chinese 
than in EA populations (0.13 versus —0.14 kg/m^ per allele) [5]. 
Even though the effect of ri?4P.?5-rs987237 on BMI was non- 
significant, the BMI-increasing allele of this SNP was marginally 
associated with increased trunk fat percentage in our population 
(P= 0.065). This might suggest that TFAP2B affects body fat 
distribution mainly through increasing trunk fat percentage. 
Although the mechanism of most genetic loci predispose 
individuals to obesity is not clear so far, evidence indicated that 
they were involved in energy uptake {FTO, MC4R, PCSKl, and 
BDMF) and adipocyte differentiation [FTO, TMFMIS, BDJVF, 
MTCH2 and NEGRI) [32-35]. 

The GRS based on 28 SNPs assesses the overall genetic 
susceptibility to obesity; each additional risk-allele was associated 
with a 0.11 kg/m^ increase in BMI (equivalent to ~320 g/allele 
for a person of 1.7 m), which tended to be smaller than the 
previously reported 0.17 kg/m^ increase of BMI per risk allele in 
~8,000 Europeans based on 32 loci [7]. Another study reported 
that each additional BMI-increasing allele in the GRS based on 
12 SNPs was associated with a 0.15 kg/m^ increase of BMI in 
~20,000 Europeans [23]. Consistently, the GRS explained less 
variance of BMI in our populations (0.90%) than in European 
ancestry populations (1.45%) [7]. Furthermore, each additional 
BMI-increasing allele increased the risk of obesity by 6% in the 
present study but approximately 10% in European ancestry 
populations [23,36]. These differences maybe due to the fact that 
the majority of SNPs (24 of 28) in the score were identified in 
European ancestry at genome-wide significant levels, whereas only 
1 1 of the SNPs did so in East Asians. Besides, smaller minor allele 
frequency of 19 SNPs in Chinese population than in Europeans 
and Winner's curse effect might be important contributions to the 
cross-population differences [37]. Nevertheless, the conclusions are 
expected since it is well established that the SNPs from GWAS 
have limited predictive ability to distinguish individuals with 
obesity from ones with normal weight [38-39]. 

We subsequentiy compared the effect of the GRS across 
adiposity and body composition traits and observed that GRS 
showed the strongest association with BMI, followed by trunk fat 
and overall body fat percentage. This may due to that the GRS is 
consisted of loci identified in GWAS for BMI rather than for body 
composition traits. Although BMI is a good marker for general 
adiposity, it does not distinguish between lean and fat body mass, 
as well as regional fat distribution. 

Contrary to our hypothesis that the EAA GRS, based on East 
Asian-identified SNPs, would have a larger effect on BMI in our 
population than the EA GRS, based on European-identified 
SNPs, the effects of EA GRS on BMI and obesity risk were largely 
similar to that of EAA GRS. Interestingly, effects of EAA GRS on 
other adiposity traits tended to be smaller than those of EA GRS. 
The reason for this discrepancy is not clear. However, we 
speculate that this might be related to the fact the BMI represents 
a somewhat different adiposity phenotype in East Asian than in 
European ancestry populations, such that BMI-associated loci 
identified in European ancestry populations are, to some extent, 
different from BMI-associated loci identified in East Asian 
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Table 4. Effect of interaction between GRSs and PA on BIVll. 





BMI 

(kg/m^) 


Main effect of GRS In 
low PA 


Main effect of GRS in 
moderate PA 


Main effect of GRS 
In high PA 


GRSxPA 
interaction 


Power 


N 


219 


1226 


1449 






GR5 


Beta 
(95% CI) 


0.20 (0.08) 


0.12 (0.03) 


0.08 (0.03) 


-0.06 (0.03) 


0.25 


P 


0.009 


0.0003 


0.004 


0.022 




EA GRS 


Beta 
(95% CI) 


0.15 (0.08) 


0.12 (0.04) 


0.09 (0.03) 


-0.05 (0.04) 


0.20 


P 


0.071 


0.0008 


0.005 


0.13 




EAA GRS 


Beta 
(95% CI) 


0.22 (0.12) 


0.15 (0.05) 


0.09 (0.04) 


-0.08 (0.05) 


0.37 


P 


0.057 


0.003 


0.038 


0.044 





Interaction beta is the difference in trait of GRS with every increase in PA category, e.g., an interaction beta of —0.06 kg/m^ for BMI represents a 0.06 kg/m^ attenuation 
in per BMI-increasing effect of GRS with every increase in PA category. Pmteraction was adjusted for age, age^, sex, region and the first two principle components. 
doi:l 0.1 371 /journal.pone.0091 442.t004 



populations. Although no heterogeneity was found between efFects 
of EA and EAA GRS on adiposity traits, this may attribute to the 
limited sample size of our current study, which provided a too 
wide confidence interval to test the significance for efFects 
heterogeneity between EAA GRS and EA GRS. In this regard, 
larger studies with detailed adiposity phenotypes will be required 
to confirm these observations. 

Consistent with observations in European ancestry populations 
[17,20,23,26], we found that the effect of GRS on obesity- 
susceptibility was more pronounced in individuals with low PA 
than in those with high PA levels. While each additional risk allele 
increases the BMI by 0.08 kg/m^ in physically active individuals, 
the increase was 0.20 kg/m^ per allele in physically inactive 
individuals, or a difference of 60%, which is twice the difference 
seen in European ancestry [17,20,23,26]. More mechanism studies 
on how PA affects the genetic association are needed. 

We acknowledge that the sample size in the present study had 
enough powers (s99%) for detecting associations of the GRSs 
with obesity and related traits, but less than 80% power for 
detecting betas, as previously reported, of most individual SNPs 
(except for PCSKl -rs26l967 and FrO-rs9939609) for associations 
with BMI, and for the gene xphysical activity interactions in 
obesity. We, therefore, cannot rule out the possibility of false 
negative results for these analyses and more study with large 
sample size is required for making firm conclusion. 

In conclusion, the effect size and explained variance of GRS on 
BMI tends to be lower in Chinese Hans than European ancestry, 
and the effects of GRSs based on all loci, EA-identified loci or 
EAA-identified loci, are the same for BMI and obesity risk. 
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